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ABSTRACT 

Hierarchical theories of structure formation predict that clusters of galaxies should be embedded in a web 
like structure, with filaments emanating from them to large distances. The amount of mass contained 
within such filaments near a cluster can be comparable to the collapsed mass of the cluster itself. Diffuse 
inf ailing material also contains a large amount of mass. Both these components can contribute to the 
cluster weak lensing signal. This "projection bias" is maximized if a filament lies close to the line-of- 
sight to a cluster. Using large-scale numerical simulations of structure formation in a A-dominated 
cold dark matter model, we show that the projected mass typically exceeds the actual mass by several 
tens of percent. This effect is significant for attempts to estimate cluster masses through weak lensing 
observations, and will affect weak lensing surveys aimed at constructing the cluster mass function. 

Subject headings: Galaxies-clusters, cosmology-theory 



1. INTRODUCTION 

Clusters of galaxies are excellent cosmological probes. 
Their size suggests that their constituents provide a fair 
sample of the universe. Their structure and hydrodynamic 
state provide information on their formation and evolu- 
tion, and thus upon models of structure formation. Mea- 
surements of the abundance of clusters of a given mass 
allows constraint of the amplitude of mass fluctuations in 
the universe; measurements of abundance evolution can be 
used to constrain the mass density f2 ro . 

Many of these approaches depend upon some knowledge 
of the mass, or mass distribution, of the cluster. Most tech- 
niques for measuring cluster masses are based upon some 
equilibrium assumption which relates the cluster mass to 
an observable such as the temperature of the intraclus- 
ter plasma or the velocity dispersion of cluster galaxies. 
Recently, however, it has become feasible to measure the 
surface density distribution of a cluster through obser- 
vations of weak gravitational lensing of the background 
galaxy field by the cluster. An attractive quality of this 
technique is that no assumptions about the dynamical or 
thcrmodynamical state of the cluster components need be 
made. In other words, weak lensing analyses probe the 
mass distribution directly. 

However, analyses of the mass distribution of a clus- 
ter drawn from weak lensing observations are n ot without 
problems (for a recent review, see Mellier 1999). Many of 
these relate to details of the procedure adopted to go from 
the observed ellipticity distribution to the mass, or from 
instrumental effects. We will not discuss these in this Let- 
ter. We are interested here in the degree to which weak 
lensing mass estimates of clusters are affected by lensing 
from material outside but nearby the cluster. This is a 
source of systematic error which is not well understood 
(thoug h it ha s been alluded to in earlier work , e.g. Miralda- 
Escude [l99lt Wambsganss, Cen & Ostriker p9g| ). 



Since clusters form in overdense regions, the volume sur- 
rounding a cluster is li kely to contain infalling overdense 
matter (Gunn & Gott 1972| ). This infalling matter could 
add to an observed lensing signal and result in an over- 
estimate of the cluster mass. It is possible that this bias 
could be quite severe. In modern hierarchical models of 
structure formation, such as the Cold Dark Matter (CDM) 
model, numerical simulations imply that clusters form pri- 
marily at the intersections of filaments in a web of cosmic 
structure, accreting additional diffuse mass and smaller 
collapsed objects that drain along these filaments. It is 
thus reasonable to expect a beaded filamentary structure 
surrounding most clusters of galaxies. Tentative observa- 
tional evidence of filamentary structure nea r clus ters has 
been report ed recently (Kull & Boehringer 1999] , Kaiser 



et al. 1999| ). A filament lying near the line-of-sight will 



also lens the background galaxies, and therefore contribute 
spuriously to the lensing signal. If the observed lensing 
signal were attributed solely to the cluster, the inferred 
cluster mass could be much larger than its actual mass. 

In this Letter, we use mock clusters from numerical sim- 
ulations to explore the significance of the systematic mass 
overestimation induced by the additional lensing signal of 
both diffuse and filamentary material near the cluster. We 
find that this effect can be quite significant and must be 
considered when evaluating the results of lensing mass re- 
construction techniques. In the next section, we describe 
the numerical data and our procedure for evaluating the 
errors introduced into cluster mass estimates by nearby 
material. Section 3 describes the results of our analyses. 
We discuss these results and outline plans for future study 
at the end. 

2. METHOD 

Weak lensing mass reconstruction techniques produce 
a map of shear or convergence. These are integrals of 
the mass along the line-of-sight times a projection kernel. 
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This kernel is quite wide in the redshift direction, scaling 
as DlDls/ Dg where Dl is the distance from the observer 
to the lens, Ds from the observer to t he sou rce and Dls 
from the lens to the source^ (Mellier 1999). Under the 
assumption that the cluster is the most massive object 
along the line-of-sight and is well localized in space (the 
thin- lens approximation), the convergence map is propor- 
tional to the projected surface density map of the lensing 
cluster itself. 

Any additional mass located near the cluster and along 
the line-of-sight will also contribute to the lensing sig- 
nal. Since the kernel is such a slowly varying function 
of distance, material even large distances from the cluster 
will contribute within the thin lens approximation. For a 
source at z = 1 and a cluster at z ~ 0.5 the kernel changes 
by only 1% within ±40/i _1 Mpc of the cluster in a uni- 
verse with f2 m — 0.3 = 1 — J1a, with similar results in 
other cosmologies. As a result, weak lensing observations 
will probe the projected density of a cluster plus all of the 
material in its vicinity. Note that this "nearby" material is 
essentially "at" the redshift of the cluster for the purposes 
of lensing, and so cannot be distinguished by using extra 
information such as source redshifts. 

To study the effect of the surrounding mass upon the 
projected mass inferred from lensing observations of the 
simulated clusters, we have examined the mass distribu- 
tion around several clusters of galaxies extracted from a 
large cosmological simulation. The simulated clusters were 
taken from the X-Ray Cluster Data Archive of the Lab- 
oratory for Computational Astrophysics of the National 
Center for Supercomputing Applications (NCSA) , and the 
Missouri Astrophysics Research Group of the University 
of Missouri. To produce these clusters, a particle-mesh N- 
body simulation incorporating adaptive mesh refinement 
was performed in a volume 256/i _1 Mpc on a side. Re- 
gions where clusters formed were identified; for each clus- 
ter, the simulation was then re-run (including a baryonic 
fluid) with finer resolution grids centered upon the cluster 
of interest. In the adaptive mesh refinement technique, 

1 For a distribution of source distances, one takes an appropriate 
integral of this expression. 
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Table 1 

Parameters for the 3 simulated clusters discussed in this Letter. 
sphere is the size of the sphere, centered on the cluster, which wc 
consider in this work, T200 is the 3D radius defined in Eq. [I] and M200 
the mass enclosed. Af>70 and -M>io are masses cutting out particles 
above thresholds of 70 and 10 times local density respectively (see 
text). Aftot is the total mass in the sphere. 



the mesh resolution dynamically improves as needed in 
high-density regions. The "final" mesh scale at the high- 
est resolution was 15.6 /i -1 kpc, allowing good resolution 
of the filamentary structure around the cluster. Inside the 
cluster, the characteristic separation between the smallest- 

1 /3 

mass particles, given by d = (4tt r% 00 /3N) ' with N the 
number of particles inside the region, was approximately 
86 kpc for all three clusters examined h ere. The code is 
described in detail in Norman & Bryan ( |1999[ ). 

The clusters used here were extracted at z — from 
simulations of a ACDM model, with parameters f2 m = 0.3, 
Vl B = 0.026, n A = 0.7, h = 0.7, and er 8 = 0.928. In this 
Letter we observe these clusters as if they were at z = 0.5. 
In future work we plan to investigate the dependence of 
these results on cosmology and on cluster redshift. 

Since the number density of rich clusters is approxi- 
mately 4>* ~ 10 -5 Mpc -3 , the typical separation between 

them is <f)~ 1/3 ~ 40/i" 1 Mpc. This is a characteristic scale 
for filaments: volumes containing a cluster and with one- 
dimensional extent ~ 40 ft, _1 Mpc should also contain much 
of the nearby filamentary structure. Three such volumes, 
containing a rich cluster (Clusters 0, 2 and 4) as well as 
satellites and filaments, were extracted from the archive. 
We selected clusters that did not appear to be mergers or 
have a large secondary mass concentration nearby. Such 
systems might be excluded observationally from studying, 
for instance, the galaxy line-of-sight velocity distribution. 
For each volume, the "center" of our cluster was deter- 
mined using a maximum-density algorithm. As the ex- 
tracted volumes were not spherical, it was possible that 
some lines of sight could contain more mass than others 
simply by geometry. To avoid such biases we restrict our 
analysis to particles that lay within the largest sphere, cen- 
tered on the cluster, which was contained entirely within 
the extracted volume. The radii of these spheres, i? sp hcrc, 
are listed in Table 0, along with other properties of the 
clusters. Note that these radii are large compared to the 
projected values of r2oo obtained for each cluster; thus no 
significant radial surface density gradient is introduced by 
a decreasing chord length through the sphere with radius. 
It is also important to note that since our volumes are by 
necessity limited, our results should be interpreted as a 
lower limit to the size of the effect; the magnitude of the 
lensing kernel is still significant at the edge of our spherical 
volume. 

Each of the clusters we examined was surrounded by a 
large amount of mass. Most of this material appeared by 
eye to be collapsed into "beads" along a filamentary struc- 
ture, although a small number of clumps could be found 
outside the filaments. We show a projection of a fraction 
of the points from the simulation of cluster 4 in Fig. [j} 
The filamentary structure and satellites are easily evident. 
Note that this filamentary structure extends well beyond 
our radius i? S phcrc- No single projection can show the full 
3D nature of the structure, in which the filamentarity is 
even more apparent. Since much of this mass is at low 
density it is unlikely it would be a site for galaxy forma- 
tion or otherwise emit light. Thus this structure would 
not be easy to constrain by observations of redshifts near 
the cluster. 
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Fig. 1. — Cluster 4, with a fraction of the particles projected onto the 
x — z plane. Note the filamentary structure, with clumps beading up in 
the filaments. The dashed circle marks the sphere of radius i? S phere to 
within which we have restricted our analysis (sec text). 

We observed each of the three selected clusters from 
10,000 randomly chosen viewing angles. For each cluster 
and viewing angle, the projected surface density map was 
constructed and used to estimate r2oo, the radius within 
which the mean interior density contrast is 200. In three 
dimensions, this radius is defined in terms of the enclosed 
mass by 

M (< r 200 ) = 200 x f il m p alt r| 00 . (1) 

The projected estimate of r2oo was extracted from the sur- 
face density map by considering the radius of the circle, 
centered on the cluster, which contained the amount of 
mass given by Eq. (|l]) above, i.e. 

/ d0 RdR E(R,0) = M(<r 200 ) (2) 
Jo Jo 

with £ (i?, 9) the surface density on the map in terms 
of a two-dimensional radius R. This radius was com- 
pared to the cluster's true r200j extracted from the three- 
dimensional mass distribution. The ratio of the projected 
mass to true mass is given simply by the cube of the ratio 
of the estimated value of r2oo to the true value. For each 
cluster, a value of this ratio was obtained for each viewing 
angle. 

3. RESULTS 

Our main result is displayed in Fig. ^, where we show the 
distribution of projected vs. "true" cluster mass in each of 
the three simulated clusters. We have checked that the 
features in the histogram do not come from shot noise 
due to discrete particles in the simulation. However, the 
"spikiness" is due to a discrete number of objects in the 
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Fig. 2. — The ratio of projected to actual mass within r2oo f° r °ur 3 
simulated clusters. In all cases nearby mass (primarily in filaments) has 
biased the projected mass distribution to higher values. The typical bias 
is a few tens of percent (b-c), with less of an effect on massive clusters 
(a). The vertical dashed line marks the mean of the distribution. 

neighborhood of the cluster. A small lump of matter near 
the cluster will project entirely within r 2 oo for a fraction 
of the lines of sight. For any such line of sight, the effect 
on the projected value of r 2 oo is identical. 

We expect the ratio M200 /M true to be greater than unity 
since only additional mass can be included in the projec- 
tion. The size of the smallest offset from unity for clusters 
and 4 is approximately twice what would be expected for 
material uniformly distributed at the mean density. This 
suggests that matter near the cluster is itself clustered and 
at higher than mean density. The width of the histogram 
in Fig. |^, as a fraction of the true mass, depends on the true 
mass of the cluster. Though we have only a few clusters 
in this study, it appears that the mass in nearby material 
is not proportional to the cluster mass, thus the relative 
effect of this material is smaller the larger the cluster. The 
total mass in the sphere, Af to t, is also listed in Table |l| for 
reference. 

The signal shown in Fig. ^ comes from (primarily fila- 
mentary) material outside the cluster and is not the well 
known projection effect arising from cluster asphericity. 
To verify this, we repeated the procedure described above 
for a subset of particles aimed at selecting the cluster 
alone. This was done by first selecting out particles with a 
local density contrast of greater than 70 (chosen because 
density profiles near r~ 2 reach a local density contrast near 
70 at a mean interior density contrast of 200); a small 
sphere containing the cluster but little nearby material 
was then cut out of this subset. The histogram produced 
by viewing the clearly prolate cluster at a large number 
of randomly chosen viewing angles produced a much nar- 
rower distribution, with a maximum offset of less than 10% 
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Fig. 3. — The profile C( r ii r 2) (see text) for cluster 4 along two 
lines of sight with extremal values of r2oo- We have used T2 = 800" 
in making this figure. Such a large radius is not currently achievable 
obscrvationally, but it minimizes the impact of objects nearby the cluster 
and provides a lower limit on the effect. Reducing this radius by a factor 
of two docs not change the result. Note that both profiles appear smooth 
and well behaved, even though they differ in mass by a factor of 1.7. 



in the mass ratio and a mean offset of approximately half 
that value. 

While it is beyond the scope of this Letter to perform a 
detailed modelling of any observational weak lensing strat- 
egy, we show in Fig. || two sample profiles obtained from 
aperture densitometry on our noiseless projected mass 
maps. Specifically, for Cluster 4, we show the profile along 
the lines of sight giving the largest and smallest r 20 o, for 
comparison. The C statistic is the mean value of the con- 
vergence, k, within a disk of radius r\ minus the mean 
value within an annulus r\ < r < r 2 (Fahlman et al. 1994, 
Kaiser 1995). Here we calculate C directly from the pro- 
jected mass, though observationally it would be computed 
from the tangential shear. We have taken r 2 = 800". Such 
a large radius is not (currently) achievable observationally, 
but it minimizes the impact of objects nearby the cluster 
and provides a lower limit on the size of the projection ef- 
fect. We have explicitly checked that reducing the radius 
to half this value does not change our result. 



In calculating the convergence «, the cluster was again 
assumed to be at a redshift of 0.5, with the lensed sources 
at a redshift of 1.0. In any real observation, of course, the 
lensed sources will span a range of redshifts. For mate- 
rial very close to the cluster, such as here, this will not 
affect our conclusions, and the error introduced by in- 
correctly modelling the redshift distribution of the back- 
ground sources is not the subject of this work. In addition 
to being easy to estimate, the ( statistic is robust and min- 
imizes contamination by foreground mass (Mellier 1999) 
because it is insensitive to the sheet mass degeneracy. This 
does not, however, remove the sensitivity to clustered ma- 
terial, as can be seen in Fig. ^. The mass which would 
be inferred from Fig. ^ along two lines of sight differ by a 
factor of 1.7. 

While the distribution of the projection effect varies 
from cluster to cluster, it seems clear that positive biases in 



projected mass of 20% are typical, with biases above 30% 
not uncommon. Furthermore, we emphasize again that 
these estimates are in fact lower limits; some lines of sight 
through the untruncated volume produced overestimates 
as large as 80% or more. While appropriate modeling of 
a mean density profile outside r 2 oo (drawn perhaps from 
simulations such as these) can be used to produce an un- 
biased estimator, the width of these histograms implies a 
large amount of scatter around such an estimator of the 
true (unprojected) mass. It is clear that this effect can be 
quite significant and must be taken into account when at- 
tempting to understand the results of mass reconstruction 
analyses. 

4. CONCLUSIONS 

Clusters of galaxies are part of the large-scale structure 
of the universe and observations of them should be consid- 
ered within this context. The filamentary structure near 
a cluster can contain a reasonable fraction of the mass of 
a cluster in tenuous material. Should a filament lie close 
to the line-of-sight to a cluster it will contribute to the 
weak lensing signal and positively bias the projected mass. 
We have shown that such a bias could easily be 30% (see 
Fig. §. 

Weak lensing remains one of the best methods for deter- 
mining the mass of clusters of galaxies. However methods 
which obtain the mass from estimates of the projected sur- 
face density must consider the effect outlined in this Let- 
ter. This is clearly of particular significance for attempts 
to determine the mass function of clusters directly through 
surveys of weak lensing-determined masses. 

We have not attempted to address the detailed question 
of how much this filamentary signal would affect a par- 
ticular reconstruction algorithm; the answer is no doubt 
algorithm dependent. We hope to return to this issue in 
future work, as well as to consider the effect of cosmologi- 
cal model and evolution with cluster redshift. 

The authors would like to thank Greg Bryan and Greg 
Daues for assistance in understanding the archive data, 
and Gordon Squires and Albert Stebbins for useful con- 
versations on lensing. This research was supported by the 
NSF. 
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